Dynamics of strain bifurcations in magnetostrictive ribbon 



o 

(N 



(N 

Q 

U 
d 



> 

O 

o 



Ritupan Sarmah and G. Ananthakrishna 

Materials Research Centre, Indian Institute of Science, Bangalore 560012, India 

We develop a coupled nonlinear oscillator model involving magnetization and strain to explain 
several experimentally observed dynamical features exhibited by forced magnetostrictive ribbon. 
Here we show that the model recovers the observed period doubling route to chaos as function of 
the dc field for a fixed ac field and quasiperiodic route to chaos as a function of the ac field, keeping 
the dc field constant. The model also predicts induced and suppressed chaos under the influence 
of an additional small-amplitude near resonant ac field. Our analysis suggests rich dynamics in 
coupled order parameter systems like magnetomartensitic and magnetoelectric materials. 
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I. INTRODUCTION 



The rich chaotic dynamics exhibited by sinusoidally 
driven nonlinear oscillators is ubiquitous to a large num- 
ber of systems such as turbulence in fluid systems [jj, 
chemical oscillators Q, and cardiac tissues Q. Manipu- 
lated properties of chaos, such as control of chaos [J] re- 
alized in driven nonlinear oscillators, are also reported in 
several disciplines, for instance, control of output of laser 
system [jj , enhancedperformance of permanent magnet 
synchronous motor [fj, and control of cardiac arrhyth- 
mias 0. Natural to multistate weakly driven oscilla- 
tors in the presence of ubiquitous noise is the stochastic 
resonance [8|, another generic noise induced cooperative 
phenomenon that enhances the signal to noise ratio. This 
phenomenon is also realized in several disciplines rang- 
ing from physics ( e.g., superconducting quantum inter- 
ference device magnetometers, optical and electronic de- 
vises) @ to biology ilQ]. However, it is rare to find such 
a broad spectrum of dynamical features in a single sys- 
tem. Surprisingly, physical realization of several of these 
features was reported by Vohra et al. flll - [l5l |. In their 
study of strain bifurcations in magnetostrictive ribbons 
subjected to the combined influence of sinusoidal (ac) 
and dc magnetic fields, the authors reported an extraor- 
dinarily large number of dynamical features such as (a) 
quasiperioidic (QP) route to chaos when the amplitude of 
the ac magnetic field h ac was increased in the presence of 
a dc magnetic field hd c , (b) period doubling (PD) route 
to chaos when \ia c was increased and the ac field was 
kept fixed [ill E3 , (c) a suppression and shift of period 
doubling bifurcation point [l2[ and induced subcritical 
bifurcation under small-amplitude near-resonant condi- 
tions [13], (d) control of chaos, specifically, suppressed 
and induced chaos with the application of a near-resonant 
perturbation to one of the subharmonics [HI, and (e) 
stochastic resonance [H]]. However, modeling such a rich 
dynamics exhibited by a single system in terms of the rel- 
evant strain and magnetic order parameters has remained 
a challenge. 

Vohra et al. emphasize that AE effect [l6j], i.e., the 
reversible change in the Young's modulus with applied 
magnetic field, is not responsible for the reported dynam- 



ics as the samples used are unannealed ribbons where the 
AE effect is insignificant, unlike the earlier reported PD 
and QP routes to chaos in the annealed metallic glass 
samples [ljj ■ Some of these features have been explained 
using appropriate normal forms [l8j by appealing to the 
universal nature of the bifurcation phenomenon. Our 
purpose is to develop a model in terms of magnetic and 
strain order parameters. In this paper we focus on three 
different experimentally observed dynamical features and 
show that the model predicts (i) the period doubling 
and (ii) quasiperiodic routes to chaos and (iii) induced 
and suppressed chaos in the presence of small-amplitude 
near-resonant perturbation. The general nature of these 
equations suggest much richer dynamics in ferromagnetic 
martensite samples that posses even stronger elastic and 
magnetic nonlinearities [19| and also in magnetoelectric 
materials [2(3]. The model also explains some old results 
on internal friction studies of martensites 211 ] . 



II. THE MODEL 

Experiments were performed using a sinusoidal mag- 
netic field h ac of frequency uj in the presence of a dc mag- 
netic field hdc- Our starting point is to write down the 
relevant free energies. In general, elastic free energy is a 
function of all components of the strain tensor. However, 
considering the fact that the samples are thin long rib- 
bons (5cm long) fixed at one end with the displacement or 
strain being monitored at the other end, it would be ad- 
equate to use a single strain order parameter e. Further, 
we work in one dimension and use dimensionless order 
parameters. We construct a minimal one-dimensional 
model with the strain order parameter e and magnetic 
order parameter m. In our model the weak magnetic 
nonlinearity is taken to drive the highly nonlinear elastic 
degrees of freedom [ll|. Then the total free energy has 
three contributions, namely, the strain free energy, mag- 
netic free energy, and magneto-elastic free energy. Thus 
the total free energy is Ft — F e + F m + F me . 

In one dimension the elastic free energy is given by 
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where the strain variable is e 
free energy density, and the second term is the grad ient 
free energy. Guided by the normal forms used flll - ll5j . 
we use a sixth-order polynomial for the elastic free energy 
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e . Here ft is a positive constant, 



but 9 can take on positive or negative values. To simplify 
further, we use A = 1. The minima of the free energy 
are located at e = ±e s with e s = [(ft+^ft 2 - (9)] 1/2 . The 
parameters ft (of the order of unity) and 9 are sufficient 
to control the minima. When 9 > 3/3 2 /4, e = is the 
true minima identified with the high-temperature phase. 
At 6 — 3ft 2 /A, we have a first order transition with the 
free energy vanishing at e 2 = ±3/3/2 and e = 0. For 9 
negative, e = state is unstable. This behavior can be 
parameterized by using 9 = = r e where Tf and T s 

are the first- and second order transition temperatures, 
respectively. 

The magnetic free energy is given by 
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where /„ 



Imag ~ T c 2 T 4 vvlu11 ' c — T c 

is the Curie temperature of the sample. The nature of 
magneto-elastic free energy is not known, particularly 
since the ribbon is a metallic glass sample with a glassy 
structure. However, coupling terms either preserve or 
break the invariance e — > — e and m —> —m in the elastic 
and magnetic free energies, respectively. In the absence 
of any information, we model it as a weighted sum of 
symmetry-preserving and -breaking terms given by 



with r c = T rp Tc , where T c 
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where £ is magnetoelastic coupling coefficient and < 
p < 1 is an adjustable weight factor. Other types of gra- 
dient coupling between e and m are ignored to keep the 
model simple. We use the Rayleigh dissipation function 
[22} to represent the damping of oscillating ribbon when 
the applied field is removed. This is given by 
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The kinetic energy is given by T = ^ J 

where u(y,t) is the displacement variable. Then, using 
the Lagrangian L = T — Ft and the Lagrange equations 

of motion 4 ( JK] - J^-r = we get 
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where L is a relevant time scale. Note that this is not 
a pure relaxational dynamics as m(y) is subject to sinu- 
soidal and dc fields. Indeed, e(y) is being driven through 
m(y), which is subject to the sinusoidal field. [Note also 
the mutual coupling between m(y) and e(y).] For a sam- 
ple that is clamped at one end, the above equations can 
be further simplified by noting that the strain at the 
free end is small even though it is in the anharmonic 
regime. Thus we assume that only the dominant mode 
of vibration is supported. We use e(y, t) — A(t) sin ky 
and m(y,t) = B{t) sinfcy, where k = tt/2L, with L rep- 
resenting the length of the ribbon. Further, we note that 
the dc field induces a finite strain value, which can be 
determined by equating the strain energy with the mag- 
netic energy. The application of the ac magnetic field 
induces oscillations around this strain value. Thus we 
use the dominant mode to represent deviation from the 
equilibrium value of the dc field, i.e., A = a — adc- Here, 
a is the strain amplitude in the presence of both mag- 
netic fields and a ( ic is that when only the dc field is 
imposed. Transforming Eqs. (|4l5l ) into rescaled variables 
(r = kt, f2 = uj/k, r' = r/fc) and using only the dominant 
mode, we get 



A(t) = 
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(7) 



The overdot now refers to the redefined time derivative. 
Equations ^ and ([7]) constitute a coupled set of nonlin- 
ear nonautonomous ordinary differential equations with 
several parameters. We shall use r e , r c , 7, T' , J7, and £ as 
free parameters. In addition, both h ac and hdc are the 
experimental drive parameters. Apart from this, we have 
no knowledge of the values of the parameters. 



III. RESULTS 

The dynamics of the above system of equations is quite 
complicated due to the presence of several parameters. 
Clearly, interesting dynamics can be expected when r e is 
in the range where three or two minima exist and when 
t c is negative. Thus we first map-out the region of the 
parameter space where interesting dynamics is seen. 

Consider the period-doubling route as a function of 



Using the equation of motion for the magnetic order hd c for a fixed h ac . Here, we keep r e = — 1 so that e = 
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FIG. 1. (Color online) Period-doubling bifurcation as a func- 
tion of hdc for r c = — 0.2, r e = — 1.0, £7 = 1,£ = 0.6,7 — 
1.592, T' = 0.09, h ac = 10.5, and p = 0.32. 




FIG. 2. (Color online) (a) Quasiperiodic orbit in the (A, A, B) 
space for h ac — 4.8 keeping r c = — 0.2, r e = — 1.0, = 
0.97, 7 = 0.0421, £ = 3.2, V = 0.09, h dc = 0.04, and p = 0.32. 

(b) The corresponding Poincare map in the (A p , A p ) plane. 

(c) Poincare map for a chaotic orbit for h ac = 7.4. 



state is unstable, although interesting dynamics is also 
observed for r e < 3(3 2 /4: as well. Other parameters are 
fixed at t c = -0.2, ft = 1, 7 = 1.592, T' = 0.09, L = 
5, £ = 0.6, and p = 0.32. 

The application of the dc field in the presence of ac 
field has two effects. First, the vibrations are centered 
around the strain induced by the dc field. This has al- 
ready been included in the definition of A. Second, the 
free energy for the magnetic order parameter is tilted to 
one side. In addition, the shape of the free energy is also 
affected due to the presence of nonlinear magnetoelastic 
coupling. In particular, note that apart from the addi- 
tive B term in Eq. ([S]), B 2 appears multiplicatively with 
A (equivalent to parametric forcing). This affects the 
shape of the effective elastic free energy. Thus we should 
expect /3 to be a function of hdc- To understand this, 
we first note that in the absence of hdc, increasing the 
amplitude of the ac field gives period-doubling bifurca- 
tion, as expected from studies on Duffing-like oscillators 
[23| . Increasing /3 also leads to period-doubling sequence 
keeping hdc = 0, but keeping h ac at a value where we find 
period one. In contrast, keeping h ac at a values where 
chaos is seen (with = 1) and increasing hd c reverses 
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FIG. 3. (Color online) Suppressed and induced chaos for S = 
10 _6 (D) and 5 — 10~ 2 (o) respectively as a function of h r . 
The inset shows the bifurcation diagram for r c = —0.2, r e = 
-1.0, £ = 1.0, n = 1,7 = 1.592, T' = 0.09, h ac = 10.5, and 
p = 0.32. 



the period doubling sequence, reflecting the compensat- 
ing relationship between (3 and hdc- This can be made 
quantitative by keeping h ac — 9.2 (a value correspond- 
ing to a period-four orbit) and increasing both hd c and (3 
such that we always see the period- four orbit. This gives 
a relation /? = 1 + 10hd c - However, while /3 is linear in 
hdc, the prefactor depends on the value of h ac used. 

For further calculations, we use this parameterized 
form of j3 in Eqs. (|6l7p . Then, increasing hd c and keeping 
h ac — 10.5 leads to period-doubling route to chaos. All 
quantitative measures are obtained after discarding the 
first 1.5 x 10 5 time steps. A plot of PD sequence is shown 
in Fig. [TJ The PD route is seen for a range of parameter 
values around those used for Fig. [1] 

Now consider quasiperiodic route to chaos. For this 
case we choose the driving frequency to be less than that 
for the period-doubling route (f2 = 1). For the results 
reported here SI = 0.97, but the similar results are ob- 
tained for a range of values of O less than unity. We 
retain most values of the parameters as in the PD se- 
quence, but use higher coupling £ = 3.2 and 7 = 0.0421. 
In this case, keeping hdc — 0.04 fixed, we sweep h ac from 
2 onwards. We find the first few harmonic frequencies (as 
found in turbulence of rotating fluid; see Ref. [l|) in the 
region h ac — 2.5 to 3.8, beyond which quasiperiodicty 
is seen until 5.1 with the emergence of two incommen- 
surate frequencies. A plot of the torus in the (^4, A, B) 
space is shown in Fig. [2a) for h ac — 4.8. The corre- 
sponding Poincare map in the (A p , A p ) plane is shown in 
Fig. [Hb). (The largest Lyapunov exponent vanishes in 
the regions of quasiperiodicity.) We find three regions of 
chaos interrupted by windows of quasiperiodicity in the 
interval h ac = 2.0 to 7.5. The first region of chaos is seen 
between h ac = 5.45 to 5.53 preceded by a phase locked 
region. This is followed by a transition to quasiperiodic 
regime (hd c — 5.53 — 5.64) followed by an abrupt tran- 
sition to chaos (hdc = 5.64 — 5.73). Another region of 
chaos is seen for h ac — 7.3 — 8.0. The Poincare map for 
the chaotic region (h ac — 7.4 ) is shown in Fig. [UJc)- 
The QP route to chaos is observed for a range of values 
of parameters around the values used for Fig. O 

We next examine the possibility of induced and sup- 
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pressed chaos when the system is subjected to small- 
amplitude near-resonant perturbation of the form h(r) = 
h r sin tl r r, with f2 r close to the first or the second sub- 
harmonic. For illustration we use fi r = Q/2—5, keeping S 
small. To do this we first locate a direct transition from 
period- two cycle to chaos as a function of hdc, keeping 
perturbing field h r = 0. We find that this transition oc- 
curs when hdc = 0.04014 (except for £ = 1.0, all other 
parameters values are the same as for PD route). A bi- 
furcation diagram for the strain amplitude is shown in 
the inset of Fig. [3] Keeping h dc — 0.04014, we apply 
the perturbing signal h r sm£l r T with 6 = 10~ 6 . Even 
for small h r we find that the onset of chaos is delayed. 
Further, the magnitude of the shift, though small, in- 
creases linearly with h r . Identifying fio with the onset of 
chaos (hd c = 0.04014) when h r = and /i with that for 
finite h r , the normalized shift (□) M ~ Mo for the onset of 
chaos (obtained by computing the Lyapunov exponent) 
as function of h r is shown in Fig. [3J We have also stud- 
ied the influence of 5 on the nature of the dynamics. For 
instance, for 8 = 10~ 2 , we find induced chaos. The mag- 
nitude of the shift in the onset of chaos (o), which is 
significantly higher than that for the suppressed chaos, 
is shown in Fig. [3] 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have developed a coupled oscillator 
model where the weak magnetic nonlinearity drives the 
highly nonlinear strain order parameter. The model ex- 
plains several experimental results on the dynamics of the 
magnetostrictive ribbon such as (a) the period-doubling 
bifurcation as a function of hdc and keeping h ac fixed, (b) 
the quasiperiodic route to chaos as a function of h ac for a 
fixed dc field, and (c) induced and suppressed chaos un- 
der the influence of the resonant perturbing field h r fill — 
fl5j . However, the suppressed chaos seen here is for small 
5 while induced chaos is seen for large S, which is the op- 
posite of what has been reported. This may be attributed 
to parametric type of forcing in our equations. Note that 
the method of control of chaos is different from tradi- 
tional methods [3] . The model demonstrates induced and 
suppressed chaos in the presence of near-resonant per- 
turbation. The model exhibits rich dynamics including 
transient chaos. For instance, for the symmetry restoring 
crisis ??, we find that the mean time < r > spent in the 
precrisis attractor scales as < t >~ (h ac — h ac )~ n , where 
h ac is the critical value, with rj ~ 1.06. 

The approach is clearly applicable to any coupled order 
parameter system fl9l |20| such as magnetomartensites 



1 19;] , which exhibits high elastic strain and magnetization, 
and ferroelectromagnetic materials, which exhibit mag- 
netization and electric polarization [20]. Such systems 
are also described by polynomial forms of free energy, as 
used in the model. Our analysis suggests rich dynamics if 
experiments are performed in a similar geometry on these 
materials. This should encourage dynamic experiments 
in these materials where experiments are traditionally 
carried out in quasistatic conditions. 

The model has been adopted to explain some old 
unexplained results [251] in the internal friction experi- 
ments on samples of nonmagnetic martensite samples of 
Cu82.9Al14.1Ni3 reported |2ll . [26j . These experiments 
have also been carried out in a similar geometry (but 
transverse drive) near martensite transformation temper- 
ature T m , where the elastic nonlinearities are significant. 
Our equations automatically describe the twinned struc- 
ture in one dimension when the magnetic free-energy con- 
tribution is dropped [27|. The model recovers all results 
[25| including the period four cycle (which the authors 
even fail to recognize; see Fig. 2 of |21(). 

These equations in their present form describe the dy- 
namics of magnetomartensites as well [l9[ . These alloys 
undergo a first-order martensitic transformation on cool- 
ing below T = T m and are also ferromagnetic below the 
Curie temperature T c . Below T m and T c they have high 
degree of elastic and magnetic nonlinearities. Indeed, the 
application of a small amount of magnetic field induces 
large strains by the rearrangement of martensite variants 
19]. Similarly, stress also affects magnetization, thereby 
suggesting strong magnetoelastic coupling. Thus a rich 
dynamics is predicted in magnetomartensites if experi- 
ments are performed in a similar geometry. It would be 
interesting to verify this prediction. 

Finally, since the application of magnetic field can in- 
duce large strains (as high as ~ 10%), magnetomarten- 
sites are better suited for actuator applications than even 
the best magnetostrictive materials (such as Terfenol- 
D, with a strain of ~ 0.2% ). Further, these materials 
should provide more accurate control of the measurement 
of dynamic strains compared to magnetostrictive mate- 
rials [3. 
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